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1.  Introduction 

V 

One  of  the  methods  of  estimating  auroral  activity  is  to 
derive  the  hemispheric  power  (HP),  which  is  the  area- 
integrated  particle  energy  deposition  into  each  high- 
latitude  hemisphere  (Fuller- Rowell  and  Evans,  1987). 
Statistical  patterns  of  electron  and  ion  energy  influx  from 
Defense  Meteorological  Satellite  Program  (DMSP)  satel¬ 
lites  were  made  as  a  function  of  Kp  by  Hardy  et  al.  (1985, 
1989),  and  patterns  of  the  total  (electron  and  ion)  energy 
input  from  National  Oceanic  and  Atmospheric  Adminis¬ 
tration  (NO/\A)  satellites  were  made  as  a  function  of  HP  by 
Fuller-Rowell  and  Evans  (1987).  We  use  the  intersatellite- 
adjusted  global  (sum  of  both  hemispheres)  electron 
hemispheric  power  (HPe)  described  by  Emery  et  al. 
(2005,  2006,  2008)  to  investigate  the  relationship 
between  solar  cycle  variations  in  the  auroral  electron 
power  {Pc)  and  the  varying  fractions  of  the  various  types 
of  structures  in  the  solar  wind  (coronal  hole  flows, 
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transients  associated  with  coronal  mass  ejections  (CME) 
at  the  Sun,  and  slow  solar  wind)  through  the  solar  cycle 
(SC).  The  aim  is  to  explain  the  large  variability  in  HPc 
exhibited  over  three  different  ~ll-year  solar  activity 
cycles  by  Emery  et  al.  (2006,  2008). 

The  highest  cross-correlations  of  HPe  are  with  the  3-h 
geomagnetic  activity  indices  of  Kp  and  Op,  where  Op  is  a 
linearization  of  the  Kp  index.  Emery  et  al.  (2008)  showed 
that  the  cross-correlations  of  daily  average  Southern 
Hemisphere  HPe  were  0.87  and  0.83  for  daily  Kp  and  Ap, 
respectively.  (The  correlations  decreased  by  ~0.17  when 
hourly  values  were  used.)  The  high  cross-correlations  with 
Kp  and  Ap  are  not  surprising  since  the  auroral  inputs,  along 
with  solar  conductance  on  the  dayside,  create  the  auroral 
current  systems  measured  as  deviations  in  the  magnetic 
field  at  ground  stations  which  lead  to  the  Kp  index.  The 
present  study  will  examine  cross-correlations  of  hourly  Pc 
with  solar  wind  parameters  such  as  the  solar  wind 
velocity  (V4w).  the  interplanetary  magnetic  field  (IMF) 
magnitude  (|B|),  the  Hz  component,  and  various  combina¬ 
tions  in  the  context  of  different  solar  wind  structures. 

The  ultimate  driver  of  the  auroral  precipitation  is  the 
Sun  via  the  solar  wind.  The  solar  wind  brings  *frozen-in’ 
magnetic  field  lines  from  the  Sun  to  the  Earth.  If  the  IMF  is 
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directed  southwards  (Bz  negative),  or  anti-parallel  to  the 
Earth’s  external  magnetic  field,  then  the  two  different 
magnetic  fields  can  merge  and  allow  direct  access  of  the 
solar  wind  to  the  Earth’s  magnetosphere  and  ionosphere. 
Southward  IMF  (negative  IMF  B^)  results  in  increased 
aurora,  stronger  polar  cap  potential  drops,  and  stronger 
currents  in  the  ionosphere  (e.g.  Reiff  and  Luhmann,  1986). 

Magnetic  reconnection  is  the  way  mass,  momentum 
and  energy  are  transferred  from  the  solar  wind  and  IMF 
into  the  magnetosphere-ionosphere-atmosphere  system. 
Most  reconnection  is  quasi-stationary  on  the  dayside  of 
the  magnetopause  (e.g.  Raeder,  2006),  but  some  recon¬ 
nection  is  time-dependant  and  patchy,  such  as  that 
experienced  by  flux  transfer  events  (FTE)  under  strong 
southward  IMF  conditions  (Bz  negative)  in  the  winter 
hemisphere  due  to  the  dipole  tilt  at  solstice  between  two 
X-lines  (Raeder,  2006).  Other  FTEs  can  develop  with 
strong  IMF  Bx  or  By  at  other  times  in  other  locations. 
Observations  by  multiple  satellites  (e.g.  Hasegawa  et  al., 
2006;  Lui  et  al.,  2008)  show  these  FTEs  as  frozen-in 
twisted  bundles  of  magnetic  fields  or  ropes  that  can  be  as 
wide  as  the  Earth,  which  last  only  a  few  minutes  as  they 
deliver  the  solar  wind  particles  flowing  along  these  ropes 
to  the  magnetosphere. 

The  cold  solar  wind  electrons  from  FTEs  or  from 
regular  quasi-steady  reconnection  on  the  dayside  magne¬ 
topause  are  then  energized  in  the  Earth’s  magnetotail  and 
are  deposited  in  the  high-latitude  auroral  regions  in 
upward  field-aligned  currents  (FACs)  (e.g.  lijima  and 
Potemra,  1976).  Discrete  electrons  are  accelerated  down¬ 
wards  within  Region  1  FACs  in  the  pre-midnight  sector 
towards  the  poleward  edge  of  the  aurora,  while  the  more 
common  diffuse  electron  aurora  caused  by  hot  plasma 
sheet  electrons  are  pitch-angle  scattered  into  the  loss  cone 
by  whistler  mode  waves  near  the  equatorial  plane  and 
precipitated  into  mostly  the  post-midnight  sector  (e.g. 
Hamrin  et  al.,  2005),  creating  the  Region  2  upward  FACs 
(e.g.  Shue  et  al.,  2001).  The  time  needed  to  energize  the 
auroral  electrons  in  the  Earth’s  magnetosphere  is  about 
1  h,  since  cross-correlations  of  the  estimated  auroral 
energy  flux  and  the  IMF  or  solar  wind  velocity  are 
maximized  with  a  1  h  lag  (see  Section  3.4). 

The  most  common  period  found  in  solar  wind  and 
geomagnetic  activity  parameters  is  ~27  days,  or  the 
equatorial  rotation  period  of  the  Sun  as  viewed  from  the 
Earth  (e.g.  Svalgaard  and  Wilcox.  1975;  Gosling  et  al.,  1976, 
1977).  Another  prominent  period  is  associated  with  the 
-•'ll  -year  solar  cycle  (SC).  Each  SC  starts  at  solar  (sunspot) 
minimum,  although  the  solar  magnetic  dipole  field 
reverses  just  after  solar  maximum.  The  change  in 
magnetic  polarity  discovered  by  Hale  and  Nicholson 
(1925)  produces  a  22-year  solar  magnetic  cycle  (the  Hale 
cycle),  which  starts  with  an  even  SC  and  ends  with  an  odd 
SC.  javaraiah  et  al.  (2005)  present  evidence  of  the 
equatorial  rotation  rate  and  the  latitudinal  gradient  of 
the  rotation  rate  from  sunspots  that  suggest  that  the  basic 
physical  cycle  of  the  Sun  is  the  22-year  magnetic  cycle 
rather  than  the  11 -year  activity  cycle.  Cliver  et  al.  (1996) 
showed  that  the  sunspot  numbers  are  highest  in  the 
ending  odd  SCs  in  the  even-odd  Hale  pairs  over  the  last 
150  years  starting  with  SC  10  in  1857.  However,  this 


pattern  was  broken  in  SC  23  as  shown  in  the  bottom  panel 
of  Fig.  la,  when  the  sunspot  number  in  SC  23  was  less 
than  that  for  SC  22  in  the  most  recent  Hale  cycle.  The  last 
time  this  happened  was  in  the  Hale  cycle  for  SC  8 
(maximum  in  1837)  and  SC  9  (maximum  in  1848), 
although  it  happened  also  in  the  4-5  Hale  cycle  and 
several  times  before  1760  when  the  sunspot  data  quality  is 
‘questionable’  or  ‘poor’  (e.g.  Cliver  et  al.,  1996). 

Our  study  investigates  the  periodicities  in  Vjw  ^nd  in  the 
global  auroral  electron  power  (Pe).  The  differences  and 
similarities  of  Pe  in  odd  and  even  SCs  are  also  examined  and 
compared  with  previous  studies  of  the  22-year  variation  in 
the  geomagnetic  activity  index  aa  by  Cliver  et  al.  (1996), 
who  showed  maximum  aa  activities  between  declining 
even  cycles  and  ascending  odd  cycles.  However,  like  the 
sunspot  numbers,  this  aa  pattern  was  broken  in  the  latest 
Hale  cycle,  when  the  yearly  aa  index  reached  a  maximum  of 
38  in  2003  during  the  declining  phase  of  SC  23.  This  yearly 
aa  maximum  in  2003  exceeded  all  other  yearly  values  in  the 
last  four  solar  cycles  (Fig.  1  of  Kozyra  et  al.,  2006). 

In  this  study,  we  consider  the  solar  wind  to  be 
composed  of  three  main  types  of  flows:  slow  speed, 
coronal  hole  (CH)  high-speed  streams  (HSS),  and  transi¬ 
ents.  The  slow-speed  wind  velocity  (Usw<'^400km/s) 
comes  from  the  borders  of  coronal  holes  (e.g.  Maris  and 
Maris,  2005)  and  from  the  region  of  the  equatorial 
streamer  belt.  Coronal  holes  are  found  at  both  high  and 
low  latitudes  on  the  Sun,  where  the  polar  CHs  are  more 
prominent  during  solar  minimum.  The  low-latitude  CHs 
give  rise  to  the  corotating  HSS  seen  at  Earth,  which  are 
particularly  prominent  in  the  descending  and  minimum 
phases  of  the  solar  cycle.  Typically,  there  are  two  CHs  of 
opposite  polarity,  which  give  rise  to  a  two-stream 
magnetic  sector  structure  (IMF  Bx  positive  and  negative) 
in  the  ecliptic  with  magnetic  activity  on  Earth  recurring 
every  27  days  and  also  every  --13.5  days  because  of  the 
two  streams.  However,  there  are  periods  when  the?c  are 
more  than  two  CHs,  such  as  described  by  Temmer  et  al. 
(2007)  where  three  CHs  produced  strong  9-day  periodi¬ 
cities  in  V^sw  in  2005. 

The  third  type  of  flows  is  associated  with  CME.  These 
include  interplanetary  coronal  mass  ejections  (ICMEs,  e.g.. 
Cane  and  Richardson,  2003  and  references  therein)  and 
the  sheath  of  compressed  solar  wind  formed  upstream  of 
many  ICMEs  that  is  often  bounded  by  a  shock  at  the 
leading  edge.  Magnetic  clouds  are  a  subset  of  ICMEs 
where  the  IMF  magnitude  is  large  and  the  field  direction 
changes  slowly  through  a  large  angle.  Magnetic  clouds  are 
often  responsible  for  the  most  intense  geomagnetic 
storms  (e.g.,  Zhang  et  al.,  2007).  Li  and  Luhmann  (2004) 
found  that  the  occurrence  frequency  of  bipolar  magnetic 
clouds  changes  from  north  to  south  in  even  SCs  to  south  to 
north  in  odd  SCs.  where  the  polarity  change  occurs  late  in 
the  descent  phase  of  the  SC. 

According  to  Berdichevsky  et  al.  (2005),  -8%  of  ICMEs 
between  1996  and  2002  were  embedded  in  HSS.  More 
commonly,  ICMEs  are  found  near  HSS  leading  edges,  or  in 
slow  interstream  solar  wind.  Watari  and  Watanabe  (2000) 
reported  about  half  of  the  magnetic  disturbances  (mini¬ 
mum  Dst<-50nT)  in  the  solar  minimum  years  of 
1995-1997  were  associated  with  ICMEs,  and  that  almost 
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3  Percent  of  Time  in  Solar  Wind  Structure 


b  Percent  of  Auroral  Electron  Power  in  Solar  Wind  Structure 


Fig.  1.  (a)  The  fraction  of  time  that  is  spent  in  each  kind  of  solar  wind  structure  from  January  1972  to  August  2008  for  hourly  average  values  plotted  as  27- 
day  averages.  Only  hours  with  both  Vsw  and  IMF  data  are  considered.  The  percentages  in  the  labels  give  the  average  fraction  of  time  spent  in  each  structure 
type  for  27-day  periods  where  the  structure  is  present,  or  for  all  times.  The  monthly  smoothed  sunspot  number  is  plotted  with  dotted  lines  to  delineate 
each  phase  of  the  solar  cycle  (minimum,  ascent,  maximum,  and  descent),  where  the  means  in  each  phase  are  plotted  as  symbols  linked  with  straight  lines, 
(b)  Same  as  (a)  for  the  fraction  of  deposited  hourly  by  specific  solar  wind  structures  from  November  1978  where  only  hours  with  HP^  and  IMF  are 
considered.  P,  lags  Vsw  and  IMF  by  1  h.  The  ‘uncertain’  type  of  solar  wind  structure  is  plotted  at  the  bottom. 


half  of  these  ICMEs  were  propagating  through  HSS.  All  of 
their  stronger  storms  (minimum  Dst<-100nT)  were 
either  definitely  or  possibly  associated  with  ICMEs. 
Richardson  et  al,  (2001)  reported  that  97%  of  the  major 


Kp  storms  (/Cp>8-)  from  1972  to  2000  were  associated 
with  ICMEs,  while  half  or  more  of  any  storms  with  Kp 
5-  or  higher  were  associated  with  ICMEs  during  solar 
maximum. 
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The  present  study  examines  the  relationship  between 
the  hourly  solar  wind  structures  and  the  electron  auroral 
precipitation  as  described  by  the  intersatellite-adjusted 
electron  hemispheric  power  (HPe)  from  Emery  et  ah  (2005, 
2006,  2008)  summed  over  both  hemispheres  (Pe)  Cross¬ 
correlations  are  found  between  various  structure  types,  Pe, 
Vsw  3nd  IMF.  The  periodicities  in  Pe  during  various  parts  of 
the  SC  are  examined  to  help  explain  the  source  of  Pe 
deposition  over  time  and  the  large  variety  seen  in  Pe  over 
three  solar  cycles. 

2.  Data  analysis 

Using  the  total  statistical  patterns  of  HP,  David  Evans 
(private  communication)  computed  HP  from  single 
satellite  passes.  This  calculation  was  adapted  to  estimate 
just  the  electron  hemispheric  power  (HPe)  from  both 
NOAA  and  DMSP  satellites  (David  Evans,  private  commu¬ 
nication  and  William  Denig,  private  communication).  Of 
course,  imagers  with  most  of  the  auroral  oval  in  their 
field-of-view  can  make  the  best  HPe  estimates  (e.g. 
Lummerzheim  et  al.,  1997;  Hubert  et  al.,  2002;  Zhang 
and  Paxton,  2008;  Luan  et  al.,  2009).  However,  the  longest 
dataset  of  HPe  estimates  is  from  single-pass  NOAA  and 
DMSP  satellites  that  have  been  intercalibrated  and 
cleaned  up  by  Emery  et  al.  (2005,  2006,  2008).  We  use 
HPe  estimates  instead  of  HP  estimates  because  only  HPe 
estimates  are  available  from  DMSP  satellites.  Ions  con¬ 
tribute  ~10%  to  HP  from  NOAA  satellites  (e.g.  Hubert  et  al., 
2002;  Emery  et  al.,  2008)  and  have  different  character¬ 
istics  from  electrons. 

We  use  average  hourly  estimates  of  HPe  from  both 
hemispheres  from  day  306  (2  November)  1978  (78306)  to 
day  225  (12  August)  2008  (08225).  Pe  is  the  sum  of  HPe  in 
the  Southern  Hemisphere  and  in  the  Northern  Hemi¬ 
sphere.  When  one  hemisphere  is  absent,  Pe  is  estimated  as 
twice  the  HPe  from  the  available  hemisphere.  Pe  is  used  in 
order  to  eliminate  the  seasonal  variations  found  in  HPe 
(Emery  et  al.,  2008),  where  the  summer  HPp  values  are 
larger  than  the  winter  values  for  low-magnitude  HPe, 
while  the  winter  HPe  values  are  larger  than  summer  for 
moderate  and  high  magnitudes  of  HPe. 

The  magnitude  of  HPe  or  Pe  used  in  the  present  study 
could  be  increased  by  a  factor  of  two  since  the  Emery  et  al. 
(2008)  HPe  estimates  as  a  function  of  Kp  are  approximately 
half  the  Zhang  and  Paxton  (2008)  TlMED/GUVl  image  HP 
estimates  in  their  Fig.  8.  Zhang  and  Paxton  (2008)  show 
that  DMSP  electron  energy  fluxes  agree  with  GUVl  images 
at  the  location  of  the  satellite  crossings,  but  that  the  DMSP 
HP  estimates  are  smaller  than  the  GUVl  estimates.  They 
discuss  several  possible  reasons  for  this,  but  another 
reason  could  be  that  DMSP  or  NOAA  averaged  energy  flux 
patterns  are  low  because  of  averaging.  However,  the 
variations  in  Pe  over  time  should  be  good,  even  if  the 
magnitudes  are  relatively  low. 

We  compare  the  Pe  data  with  a  classification  of  the 
hourly  solar  wind  structures  at  Earth  during  1972-2008, 
which  is  an  update  of  that  used  by  Richardson  et  al.  (2000, 
2002).  Near-Earth,  1-h  averaged  solar  wind  plasma  and 
magnetic  field  data  from  the  Space  Physics  Data  Facility 


OMNI  database  are  used  to  classify  the  solar  wind  into 
three  types  of  structures; 

( 1 )  Transient  structures,  in  particular  ICMEs,  believed  to  be 
the  counterparts  of  CMEs  at  the  Sun  and  the  associated 
upstream  sheaths  of  compressed  solar  wind.  The  in-situ 
signatures  of  ICMEs,  which  can  help  to  identify  such 
structures  are  summarized,  for  example,  by  Zurbuchen 
and  Richardson  (2006).  For  cycle  23,  we  base  the 
classification  on  an  updated  and  extended  version  of 
the  ICME  list  of  Cane  and  Richardson  (2003),  but  also 
include  sheath  flows  following  interplanetary  shocks 
that  do  not  include  an  encounter  with  an  ICME. 
Interplanetary  shocks  are  more  extended  than  the 
driving  ICMEs  (e.g.,  Borrini  et  al.,  1982;  Cane,  1988). 
For  earlier  cycles,  similar  identifications  have  been 
made  and  used  in  various  studies  (e.g..  Cane  et  al.,  1996; 
Richardson  and  Cane,  1993,  1995;  Richardson  et  al., 
1997).  Jian  et  al.  (2006b)  have  also  compiled  lists  of 
ICMEs  and  discussed  their  identification  criteria. 

(2)  Corotating  high-speed  solar  wind  streams  from 
coronal  holes.  As  described  by  Belcher  and  Davis 
(1971),  these  HSS  are  characterized  by  an  increase  in 
solar  wind  speed  at  the  leading  edge  over  the  order  of 
a  day,  followed  by  a  slower  decay  lasting  for  one  to 
several  days.  The  interaction  of  the  stream  leading 
edge  with  the  slower  preceding  solar  wind  generates  a 
region  of  dense  plasma  with  enhanced  magnetic 
fields,  the  corotating  interaction  region  (CIR),  at  the 
stream  leading  edge. 

(3)  Interstream,  slow-speed  (<--400km/s)  solar  wind. 

A  fourth  category  includes  periods  in  which  the 
classification  is  ’uncertain’  because  of  insufficient  or 
unavailable  solar  wind  observations.  This  is  a  particular 
problem  between  late  1982  and  1994  when  the  OMNI 
data  are  predominantly  from  IMP  8,  which  was  located  in 
the  solar  wind  for  only  ^60%  of  each  orbit. 

On  occasion,  other  data  sources  can  fill  in.  For  example, 
geomagnetic  sudden  storm  commencements  (SSC)  can 
help  identify  interplanetary  shocks  that  are  often  asso¬ 
ciated  with  ICMEs,  while  cosmic  ray  depressions  can 
indicate  the  presence  of  ICMEs  and  HSS  (e.g.,  Richardson 
et  al.,  1996).  For  further  details  of  the  classification,  see 
Richardson  et  al.  (2000).  In  addition  to  the  studies  noted 
above,  Jian  et  al.  (2006a)  have  catalogued  what  they 
describe  as  “Stream  interaction  regions’’  in  solar  cycle  23. 
For  earlier  cycles.  HSS  are  listed  by  Lindblad  et  al.  (1989) 
and  Mavromichalaki  and  Vassilaki  (1998),  but  these 
include  both  corotating  HSS  from  coronal  holes 
(CH-HSS)  and  transient  flows,  the  latter  listed  as  ’’flare- 
generated”  high-speed  streams  (FG-HSS). 


3.  Solar  wind  structures  and  auroral  electron  power 

3.1.  Contribution  of  different  types  of  solar  wind  structure  to 
auroral  precipitation 

Fig.  la  shows  the  time  spent  in  each  solar  wind 
structure  using  hourly  solar  wind  structure  classifications 
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from  1  January  1972  to  2  August  2008  for  hours  when  Vsw 
and  IMF  are  available  (276611  h  or  69%  of  the  time)  plotted 
as  27-day  averages.  The  two  percentages  written  above 
each  graph  are  first  the  mean  of  the  plotted  27-day 
averages,  while  the  second  percentage  represents  the 
fraction  of  time  that  type  of  solar  wind  structure  is 
present  for  all  available  27-day  periods,  where  ‘available’ 
excludes  27-day  periods  of  missing  data,  but  includes 
periods  when  the  specified  structure  is  absent.  The  biggest 
difference  in  percentages  is  for  the  transient  flows  because 
transients  are  absent  in  29%  of  the  available  27-day 
periods,  thus  dropping  the  total  fraction  of  time  transients 
are  present  from  24%  to  17%  of  all  times. 

The  bottom  plot  shows  the  smoothed  monthly  sunspot 
number.  We  use  the  definitions  of  Maris  and  Maris  (2005) 
for  SCs  20-22,  where  the  minimum  (N)  is  defined  by  the 
months  where  the  monthly  smoothed  sunspot  number  is 
less  than  20  and  the  ‘shape’  of  the  maxima  (X)  is  included 
above  varying  limits  of  the  monthly  smoothed  sunspot 
number.  Thus,  we  define  the  SC  23  maximum  to  be 
encompassed  in  the  interval  with  monthly  smoothed 
sunspot  numbers  greater  than  103.  The  ascent  (A)  and 
descent  (D)  phases  of  the  SC  are  the  intervals  in-between 
the  minima  and  maxima.  These  four  phases  of  each 
sunspot  cycle  are  denoted  by  dashed  lines  in  Fig.  la.  Note 
that  the  plots  in  this  study  are  of  27-day  averages  of 
hourly  averages,  but  because  each  year  has  13.5  27-day 
periods,  the  14th  ‘27-day’  average  in  a  year  is  actually  14 
or  15  days  instead  of  27  days. 

We  find  averages  (excluding  ‘available’  27-day  averages 
when  the  specified  structure  was  absent)  of  the  27-day 
averages  that  are  approximately  within  the  monthly 
periods  for  each  part  of  the  solar  cycle  using  begin  and 
end  year  and  day  numbers  in  Table  1.  These  SC  phase 
averages  are  plotted  as  symbols  with  connecting  lines  for 
transients,  HSS  and  slow-speed  regimes  in  Fig.  la.  The 
uncertain  category  is  not  plotted  in  Fig.  la,  but  is  shown 
in  Fig.  lb. 

As  previously  noted,  (e.g.,  by  Richardson  et  al.,  2000, 
2002)  the  fraction  of  the  solar  wind  occupied  by 
transients  follows  the  solar  activity  cycle.  The  average 
transient  fraction  of  27-day  periods  with  transients 
present  is  24%,  whereas  the  average  for  all  ‘available’  27- 
day  periods  is  17%.  Transients  are  present  for  ~37%  of  the 
‘available’  hours  in  the  maxima  of  SC  21  and  SC  22, 
dropping  to  ~31  %  in  SC  23.  At  solar  minimum,  the  average 
occurrence  frequency  of  transients  is  -^5%. 


HSS  from  coronal  holes  are  almost  always  present 
during  each  27-day  interval,  and  are  only  absent  (as 
opposed  to  ‘missing’  where  no  data  are  available)  ~4%  of 
the  time.  The  largest  fractions  of  HSS  are  ~52%  in  the 
descent  phase  of  the  SC  with  the  smallest  fractions  of 
~32%  in  SC  21  and  SC  23  solar  maxima  (-^24%  in  SC  22 
maximum).  The  slow-speed  solar  wind  makes  up  the 
remaining  solar  wind  structure,  with  the  largest  contribu¬ 
tion  in  solar  minima  of  *^43%  (SC  20/21  and  SC  23/24) 
to  -^48%  (SC  21/22  and  SC  22/23),  and  the  smallest 
contributions  in  solar  maxima  and  the  descending  phases 
of -30%  in  SCs  20  and  21  and  -35%  in  SCs  22  and  23. 

These  occurrence  rates  are  similar  to  those  found  for 
solar  cycles  20-23  in  Richardson  et  al.  (2002),  except  our 
solar  cycle  variation  in  the  slow-speed  wind  is  larger.  The 
slow-speed  wind  and  HSS  are  out  of  phase  with  each 
other  in  Fig.  la  because  one  is  seen  at  the  expense  of 
the  other. 

Fig.  lb  is  similar  to  Fig.  la  except  it  shows  the  percent 
of  hourly  average  auroral  electron  power  (Pe)  input  that 
can  be  attributed  to  the  different  types  of  solar  wind 
structures  from  2  November  1978  to  2  August  2008  for 
hours  when  Pe,  Vsw.  and  IMF  are  available  (183825  h  or 
71%  of  the  time).  A  1  h  lag  time  between  the  IMF  or 
and  Pe  is  incorporated  because  cross-correlations  between 
Pe  and  Vsw,  Bz,  |B|,  and  various  combinations  are  largest 
with  a  1  h  lag  time  between  the  IMF  and  and  the 
resulting  deposition  of  Pe  as  described  in  Section  3.4.  We 
summed  the  average  hourly  Pe  for  each  type  of  structure 
within  a  27-day  period,  and  divided  by  the  total  summed 
Pe  during  this  27-day  period  (including  ‘uncertain’ 
structure)  to  derive  the  27-day  fractions  that  are  plotted. 
The  SC  characteristics  are  similar  to  those  in  Fig.  la,  but 
the  magnitudes  are  different.  Thus,  the  contribution  of 
transients  over  the  entire  interval  increases  from  17%  in 
Fig.  la  to  23%  in  Fig.  lb  (—35%  increase),  and  that  of  HSS 
flows  increases  from  43%  to  47%  (—10%  increase),  while 
the  contribution  of  the  slow  solar  wind  decreases  from 
37%  to  27%  (—30%  decrease).  The  increase  in  Pc  due  to 
transients  is  because  transients  are  more  ‘geoeffective’,  as 
explained  in  Section  3.2. 

We  find  HSS  contribute  —57%  to  Pe  in  the  descent  and 
minimum  phases,  decreasing  to  —32%  in  solar  maximum. 
Transients  account  for  —43%  of  Pe  at  solar  maximum  for 
SC  21  and  SC  22,  and  —38%  for  SC  23,  with  the  average 
transient  contribution  dropping  to  —6%  in  solar  minimum. 
Slow  solar  wind  contributes  most  to  Pc  in  solar  minima 


Table  1 

Parts  of  solar  cycles  20-23. 


sc#/# 

Minimum 

‘27-Day’  min. 

Maximum 

>SS# 

*27“Day’  max. 

19/20 

January  1964-October  196S 

64001 -6S270 

December  1967-August  1970 

>98 

6732S-70243 

20/21 

April  197S-February  1977 

75082-770S4 

March  1979- December  1981 

>138 

790SS-8136S 

21/22 

February  198S-February  1987 

8S028-870S4 

December  1988-November  1991 

>138 

8832S-91324 

22/23 

May  199S-May  1997 

9S109-97162 

October  1999-june  2002 

>103 

99271-02189 

23/24 

February  2006- 

06028- 

The  solar  cycle  (SC)  changes  in  solar  minimum  with  the  rise  in  the  sunspot  number,  so  solar  minimum  periods  with  smoothed  monthly  values  less  than  20 
are  shared  between  cycles.  August  2008,  the  end  of  our  study,  is  still  in  SC  23.  We  use  the  definitions  of  SCs  20-22  from  Maris  and  Maris  (200S),  and  we 
define  solar  maximum  for  SC  23  to  include  monthly  smoothed  sunspot  numbers  >  103.  The  annual  27-day  begin  and  end  year  and  daynumber  values 
used  for  each  minimum  and  maximum  period  are  also  listed. 
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27-Day  Averages  o  Hourly  Average  Values 


Fig.  2.  (a)  The  average  hourly  solar  wind  speed  ( Vsw)  attributed  to  the  three  types  of  solar  wind  structure  plotted  as  27 -day  averages  from  January  1972  to 
August  2008  for  all  days  with  Vsw  and  IMF,  where  the  ‘uncertain’  type  is  not  plotted.  The  black  line  is  the  total  (including  ‘uncertain’  but  not  including 
‘missing’)  Vsw  The  average  wind  speeds  are  given  for  each  solar  wind  structure,  as  well  as  the  percents  of  the  total  Vjw  with  27-day  periods  where  the 
structure  exists.  The  cross-correlation  coefficient  is  between  each  structure  type  Vsw  and  the  total  as  27 -day  periods  when  the  structure  exists.  The 
structure  means  in  each  solar  cycle  phase  are  plotted  as  symbols  linked  with  straight  lines,  (b)  Same  as  (a)  for  the  average  hourly  scalar  IMF  magnitude  |B|, 
where  the  black  line  is  the  total  IMF  magnitude,  (c)  Same  as  (a)  for  the  average  hourly  incoming  auroral  electron  power  (P,)  from  November  1978  where 
only  hours  with  Pe.  V^sw.  and  IMF  are  considered,  where  Pe  lags  Vsw  and  IMF  by  1  h.  The  black  line  is  the  total  P^  including  those  hours  when  or  IMF  are 
missing  as  in  1982^1983. 
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27-Day  Averages  of  Hourly  Average  Values 
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Fig.  2.  (Continued) 


(from  -27%  in  SC  23/24  to  -38%  in  SC  22/23),  and  -24%  in 
solar  maxima  and  the  descending  phases. 

These  contributions  to  Pe  are  similar  to  those  found  by 
Richardson  et  al.  (2002)  for  the  3-h  aa  magnetic  index  in 
SCs  20-23.  They  found  HSS  contributed  —67%  to  aa  in 
solar  minimum  and  —30%  at  solar  maximum,  while 
transients  accounted  for  —50%  of  aa  at  solar  maximum 
and  —10%  in  solar  minimum.  Slow  solar  wind  contributed 
—20%  to  aa  throughout  the  solar  cycle,  with  —35% 
contribution  in  the  SC  22/23  minimum.  Mavromichalaki 
et  al.  (1988)  also  studied  HSS  and  transients  (as  FG*HSS) 
with  respect  to  the  daily  magnetic  Ap  index  in  terms  of 
cross-correlations  between  the  maximum  solar  wind 
velocity  speed  and  Ap. 

3.2.  Solar  cycle  variations  due  to  solar  wind  structure 

The  reason  for  the  increased  importance  of  transients 
and  HSS  to  Pe  compared  with  their  occurrence  rates  is 
because  they  are  associated  with  higher  Vsw.  and  with 
higher  IMF  magnitudes  than  the  slow-speed  wind. 
Brautigam  et  al.  (1991)  binned  DMSP  electron  and  ion 
energy  fluxes  as  a  function  of  magnetic  latitude  and  local 
time  for  various  Bz  and  Vsw  conditions.  They  found 
minimum  values  of  HPe  for  Bz  — +4nT.  HPe  increased 
strongly  as  a  function  of  negative  Bz,  with  only  modest 
increases  in  HPe  for  larger  positive  Bz*  HPe  also  doubled  or 
tripled  for  Ww  between  346  and  572km/s  for  the  same 
negative  or  positive  Bz  values.  These  findings  were 
repeated  by  Emery  et  al.  (2008)  using  the  single-pass 
NOAA  and  DMSP  satellite  estimates  of  HPe  instead  of 


binned  patterns.  Liou  et  al.  (1998)  studied  various  MLT 
segments  of  HPe  from  POLAR  UV  images  and  found  IMF  By, 
as  well  as  Bz  and  Vsw  variations  in  HPe. 

Figs.  2a  and  b  plot  the  Vsw  and  IMF  |B|  magnitudes  that 
are  associated  with  the  three  solar  wind  structures  from 
1972  to  2008.  The  black  lines  plot  the  27-day  averages  of 
the  average  hourly  Vsw  and  |B|,  while  the  red,  green  and 
blue  lines  plot  the  same  for  each  solar  wind  structure.  The 
average  solar  wind  speed  is  highest  for  HSS  at  499km/s, 
and  lowest  for  slow-speed  winds  at  368  km/s.  The  average 
Vsw  for  transients  is  in-between  at  450  km/s.  The  solar 
wind  speed  for  the  ‘uncertain’  category  is  not  plotted,  but 
is  part  of  the  total  V^w  in  the  black  lines  of  Fig.  2a.  The 
average  total  Vsw  values  are  largest  in  the  descending 
phase  of  the  solar  cycles  in  Fig.  2a  because  the  maximum 
occurrence  frequencies  of  HSS  (with  the  largest  Vsw 
values)  are  usually  in  the  descending  phases. 

The  cross-correlations  are  between  the  27-day 
averages  of  the  total  Vsw  ^nd  Vsw  from  different  solar 
structure  types.  The  largest  cross-correlation  (r  =  0.82)  is 
for  HSS  flows  from  coronal  holes,  although  HSS  only 
account  for  -47%  of  the  total  Vsw  magnitude.  Thus 
solar  wind  peaks  are  almost  always  defined  by  HSS 
structures.  The  HSS  structures  are  particularly  promi¬ 
nent  in  Fig.  2a  in  the  total  velocity  peaks  around  the 
descending  phases  of  all  solar  cycles.  Richardson  et  al. 
(2002)  also  concluded  that  the  mean  solar  wind  speed 
most  closely  follows  variations  in  the  speed  of  corotating 
streams  (HSS).  As  they  point  out,  the  speeds  of  CME- 
related  structures  tend  to  converge  toward  that  of  the 
background  solar  wind  (slow  speed  or  HSS)  (e.g.  Lindsay 
et  al.,  1999). 
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We  used  the  scalar  IMF  magnitude  (‘Bscalar’  or  |B|), 
rather  than  the  vector  average  (‘Bvector’),  since  the  scalar 
average  from  the  NASA  OMNI  hourly  data  set  is  composed 
of  the  average  of  all  the  IMF  magnitudes  in  a  particular 
hour.  We  found  ‘Bscalar’  to  be  ~12%  larger  than  the  hourly 
vector  average  ‘Bvector’,  which  is  the  square  root  of  the 
sum  of  the  squares  of  the  hourly  IMF  component  averages. 
Charan  Dwivedi  et  al.  (2009)  found  ‘Bscalar’  cross¬ 
correlations  with  ‘Bvector’  to  be  greater  than  0.90  on  an 
hourly  basis,  dropping  to  0.67  on  a  daily  basis  where  the 
daily  Bx.  ^nd  B^  values  are  used  to  find  ‘Bvector’, 

The  IMF  scalar  magnitude  |B|  in  Fig.  2b  shows  an 
increase  during  each  solar  maximum.  All  the  maxima 
show  the  bite-out  in  |B|  around  the  beginning  of  the 
reversal  of  the  magnetic  field  polarity  ('^lOSO,  ~1990,  and 
~2001).  which  is  related  to  the  reduction  in  the  open 
magnetic  flux  associated  with  the  solar  magnetic  field 
polarity  reversal  (e.g.  Richardson  et  al.,  2000,  2002).  This 
feature  is  relatively  short  for  SCs  21  and  23  (e.g.  Maris  and 
Maris  (2005)  give  the  periods  of  the  polarity  reversals), 
while  even  solar  cycles  undergo  longer  threefold  polarity 
reversals  (Makarov  and  Sivaraman,  1989).  The  peak 
average  values  of  the  total  |B|  occur  during  the  maximum 
and  descending  phases  of  the  SC. 

Transients  carry  much  larger  average  IMF  (9.02  nT) 
than  the  total  average  |B|  because  CMEs  often  have 
intrinsically  strong  internal  fields  (especially  magnetic 
clouds).  In  addition,  compression  in  the  upstream  sheath 
also  increases  the  field  strength  associated  with  transient 
structures.  The  hourly  |B|  of  transients  does  not  correlate 
well  (r  =  0.40)  with  the  total  |B|  when  using  27-day 
averages  of  the  hourly  values.  In  agreement  with  the 
observations  of  Richardson  et  al.  (2002),  |B|  for  both  HSS 
and  slow-speed  solar  wind  structures  correlate  well  with 
the  total  |B|,  although  the  magnitude  is  ~15%  less  in  slow- 
speed  winds. 

Fig.  2c  is  a  plot  of  Pe-  Like  Vsw.  the  total  Pe  values  are 
largest  in  the  descending  phases  of  the  solar  cycle.  Also 
like  Vsw.  Pe  peaks  are  determined  mostly  by  HSS  periods 
where  the  cross-correlation  with  the  total  average  27-day 
Pe  was  0.82,  decreasing  to  0.79  with  the  full  total  Pe 
including  periods  when  the  solar  structure  was  missing. 
HSS  sources  drive  the  magnitude  and  time  variation  of  the 
total  Pe  because  HSS  are  present  43%  of  the  time  (from 
Fig.  la),  and  because  HSS  sources  result  in  an  average  Pe  of 
33.0  GW  versus  19.8  GW  from  slow-speed  sources  (from 
Fig.  2c).  Transients  are  very  important  during  the  17%  of 
the  time  they  are  present  (from  Fig.  la),  and  they  make  up 
23%  of  the  total  Pe  from  Fig.  2c.  Transients  also  determine 
most  of  the  storm-time  variations  because  of  their  large 
|B|  (9.02  nT  from  Fig.  2b),  but  they  are  absent  83%  of  the 
time.  Therefore,  HSS  sources  determine  the  average  total 
magnitude  of  Pe,  even  though  they  represent  a  little  less 
than  half  the  solar  wind  source  types. 

There  is  a  strong  similarity  of  peaks  in  Vsw  and  peaks  in 
Pe.  The  27-day  averages  give  cross-correlations  of  0.58  for 
HSS  Vsw  and  Pe,  0.51  for  slow-speed  Vsw  and  Pe,  and  0.43 
(not  significant)  for  transient  and  Pe.  The  lack  of 
correlation  of  Vsw  and  Pe  for  transients  is  probably  due  to 
the  fact  that  the  IMF  B^  is  more  important  than  Usw  for 
transients  in  the  production  of  Pe.  Emery  et  al.  (2008) 


showed  that  HPe  varied  most  with  Bz  negative,  and 
secondarily  HPe  varied  with  Vsw  magnitude.  Section  3.4 
will  show  hourly  cross-correlations  between  various  Pe 
structures  and  solar  wind  parameters  under  negative  and 
positive  Bz  conditions. 

The  characteristics  of  the  average  Usw  and  average  Pe 
are  very  different  in  each  solar  cycle.  Maris  and  Maris 
(2005)  studied  the  ‘importance’  parameter  for  transients 
(their  FG-HSS)  and  HSS  (their  CH-HSS)  between  1964  and 
1996  (SC  20-22).  The  ‘importance’  parameter  is  defined  as 
the  duration  of  the  HSS  in  days  times  the  velocity  gradient 
from  the  peak  Vsw  during  the  HSS  and  Vsw  before  and  after 
the  HSS.  They  found  transients  to  be  most  important  in 
solar  maxima,  with  the  largest  values  in  SC  20.  The 
importance  parameters  for  transients  in  SC  21  and  SC  22 
were  similar  (as  we  see  also  in  Fig.  la),  with  a  series  of 
transient  peaks  followed  by  minima  instead  of  a  single 
increase  and  decrease  in  importance  as  in  SC  20. 

The  HSS  importance  parameter  was  roughly  four  times 
more  important  than  the  transient  importance  parameter 
over  their  study  period  based  on  annual  values  in  their 
Fig.  1.  HSS  importance  (their  CH-HSS)  was  largest  in  the 
descending  part  of  the  SC,  and  was  least  in  the  maxima. 
From  Fig.  la,  we  also  see  HSS  are  most  important  in  the 
descending  phases  and  least  important  in  the  maxima  in 
terms  of  the  occurrence  frequency.  Maris  and  Maris 
(2005)  found  the  largest  HSS  importance  values  between 
SC  20/21  minimum  and  the  descent  phase  of  SC  21,  From 
Figs  1  a  and  2a,  we  find  the  occurrence  frequency  of  HSS 
and  the  magnitude  of  Vs^  to  be  about  the  same  for  SC  21 
as  for  the  other  cycles.  Maris  and  Maris  (2005)  found  HSS 
in  SC  20  and  SC  22  to  be  similar  in  the  ascent  and  maxima 
phases.  They  suggested  with  Mavromichalaki  et  al.  (1999) 
that  the  22-year  magnetic  cycle  imprints  different 
features  in  the  occurrence  rate  of  the  HSS  in  the  ascendant 
and  descendant  phases  of  the  odd  and  even  cycles.  We 
will  discuss  the  22-year  magnetic  cycle  in  terms  of  Vsw 
from  Figs  la  and  2a. 

3.3.  22-Year  Hale  cycle  variations  in  solar  wind  and  auroral 
electron  power 

Cliver  et  al.  (1996)  showed  that  the  22-year  cycle  in 
geomagnetic  aa  activity  is  characterized  by  high  activity 
during  the  second  half  (descending  and  minimum  phases) 
of  even-numbered  solar  cycles  and  the  first  half  (ascend¬ 
ing  and  maximum  phases)  of  odd-numbered  cycles,  or  aa 
is  strongest  in  the  middle  of  each  22-year  Hale  cycle.  We 
anticipate  that  this  increased  magnitude  may  also  be 
found  in  Pe  because  of  the  close  correspondence  of  Pe  with 
3-h  Kp  and  ap  (Emery  et  al.,  2008).  Russell  and  McPherron 
(1973)  argued  that  because  of  the  7"’  tilt  of  the  solar 
equator  to  the  ecliptic  as  manifested  in  the  Rosenberg  and 
Coleman  (1969)  effect,  the  excess  positive  GSE  By  in  the 
fall  and  excess  negative  GSE  By  in  the  spring  results  in  an 
excess  negative  GSM  Bz  component  that  produces  more 
magnetic  activity  on  Earth  in  the  middle  of  each  Hale 
cycle.  Cliver  et  al.  (1996)  acknowledged  the  combined 
Russell-McPherron  and  Rosenberg-Coleman  geometrical 
effects,  but  they  reminded  the  community  that  there  is  an 


BA.  Emery  et  ai  /  journal  of  Atmospheric  and  Solar-Terrestrial  Physics  71  (2009)  1157-1175 


1165 


additional  intrinsic  solar  Hale  cycle.  The  Hale  cycle  shows 
higher  sunspot  numbers  in  the  second  half  of  even- 
numbered  SCs  and  in  the  first  half  of  odd-numbered  SCs, 
where  the  maximum  sunspot  numbers  are  found  in  the 
odd  SC  maximum.  They  proposed  a  relative  excess  of  HSS 
in  the  descent  and  minimum  of  even  SCs  in  conjuction 
with  stronger  solar  poloidal  magnetic  fields,  and  a  relative 
excess  of  transients  during  the  rise  and  maximum  of  odd 
SCs.  Since  the  Hale  cycle  was  reversed  in  SCs  22-23, 
where  the  highest  sunspot  numbers  were  in  the  even  SC 

22  instead  of  in  the  odd  SC  23,  we  might  expect  to  find  a 
relative  excess  of  HSS  in  the  descent  and  minimum  of  SC 

23  instead  of  in  the  descent  and  minimum  of  SC  22. 

From  Fig.  la,  we  see  that  the  occurrence  frequencies  of 

HSS  in  the  descending  and  minimum  phases  of  SC  20  are 
larger  than  in  SC  21,  in  agreement  with  the  hypothesis  of 
Cliver  et  al.  (1996),  although  the  HSS  and  total  V5W 
magnitudes  are  similar.  For  the  unusual  Hale  22-23  cycle, 
the  HSS  V4w  occurrence  frequencies  and  magnitudes  in 
D22  and  D23  are  similar,  but  the  HSS  occurrence 
frequency  in  the  incomplete  N23/24  phase  is  slightly 
larger  than  in  D23,  and  is  8%  larger  than  in  N22/23.  In 
addition,  the  HSS  and  total  V^w  magnitudes  for  N23/24  are 
11%  and  8%  larger,  respectively,  than  in  N22/23.  In  Fig.  lb, 
the  fraction  of  Pe  deposited  by  HSS  is  exceptionally  high  in 
the  descending  phase  of  SC  22,  but  the  fraction  of  Pe 
deposited  by  HSS  is  highest  in  the  minimum  phase  of  SC 
23/24  compared  with  the  22/23  minimum.  The  magni¬ 
tudes  of  HSS  and  total  Pe  are  '^15%  higher  in  D23 
compared  to  D22.  but  are  approximately  the  same  in 
N22/23  and  N23/24.  From  Fig.  2b.  the  total  IMF  magnitude 
is  about  the  same  in  D22  and  D23,  but  is  '-14%  smaller  in 
N23/24  compared  to  N22/23.  Since  the  magnitude  of  Pe  is 
influenced  by  IMF  as  well  as  Vsw.  the  decrease  in  IMF 
balances  the  increase  in  in  N23/24  compared  to 
N22/23  in  the  resulting  magnitude  of  Pe.  The  increased 
occurrence  frequencies  of  HSS  in  Pe  and  Vsw  in  N23/24 
compared  with  N22/23  is  opposite  to  the  normal  Hale 
cycle,  but  fits  with  the  abnormal  22-23  Hale  cycle. 

When  we  examine  the  ascending  parts  of  the  solar 
cycle,  we  see  from  Fig.  1  a  that  transients  account  for  '-25% 
of  the  solar  wind  structures  in  SCs  21  and  23,  but  only 
'-11%  of  all  times  in  SC  22,  in  agreement  with  the 
hypothesis  of  a  relative  excess  of  transients  in  the 
ascending  phases  of  odd  SCs.  This  relative  excess  of 
transients  is  also  seen  in  the  fraction  of  Pe  deposited  by 
transients  in  Fig.  1  b.  However,  the  occurrence  frequency  of 
transients  is  largest  in  the  even  SC  22  maximum 
compared  with  the  odd  SCs  21  and  23,  which  is  consistent 
with  an  abnormal  22-23  Hale  cycle.  The  results  of  the  last 
three  solar  cycles  are  mixed  with  respect  to  regular  Hale 
cycles  in  Pe  due  to  changing  ratios  of  HSS  and  transients, 
partly  because  of  the  unusual  22-23  Hale  cycle. 


3.4.  Cross-correlation  of  Vsw,  IMF,  and  Pe  for  different  solar 
wind  structures 

Similar  to  the  cross-correlations  found  for  the  southern 
hemisphere  HPe  by  Emery  et  al.  (2008),  the  highest  cross¬ 
correlations  for  hourly  Pe  between  1978  and  2008  are  with 


3-h  Kp  and  Op  at  0.72  and  0.68,  respectively,  (0.90  and  0.84 
for  daily  Pe).  Thus,  we  anticipate  that  cross-correlations  of 
Pe  with  Vsw  and  IMF  will  be  similar  to  cross-correlations 
with  daily  Ap.  Charan  Dwivedi  et  al.  (2009)  have  reviewed 
the  literature  of  cross-correlations  of  interplanetary 
plasma  and  fields  with  daily  or  longer  Ap  values.  Crooker 
et  al.  (1977)  found  6-month  and  yearly  averages  of>\p  to  be 
highly  correlated  with  (r=  0.86),  Vsw^s  (r  =  0.8),  or 
Vj^Bs  (r  =  0.9),  where  Bs  is  the  southward  component  of 
Bz,  or  Bz  negative.  Ahluwalia  (2000)  focused  on  the  IMF 
magnitude  \B\  instead  of  Bs.  ^nd  found  good  correlations 
(r  =  0.92)  with  Vsw|P|  and  Vjy^\B\  for  annual  averages. 
Charan  Dwivedi  et  al.  (2009)  found  cross-correlations  of 
daily  Ap  with  Vsw|P|  to  vary  between  0.61  in  1985  to  0.86 
in  1986  and  2001  from  their  Fig.  4,  with  relatively  low 
cross-correlations  for  even  SC  ascending  phases  and 
relatively  high  cross-correlations  for  odd  SC  ascending 
phases. 

Recently,  Lei  et  al.  (2008b)  found  11 -day  running 
means  of  the  neutral  thermospheric  density  at  400  km 
to  correlate  with  Kp  with  cross-correlations  increasing 
from  0.78  to  0.96  between  2002  and  2007  from  their 
Fig.  2.  In  the  same  figure,  they  found  cross-correlations 
with  log(V5^|B|)  to  increase  from  0.71  to  0.92  in  the  same 
time  period.  The  neutral  density  is  affected  by  Joule 
heating  in  high  latitudes  caused  by  increased  conductance 
from  particle  precipitation  (Pe)  and  by  electric  fields 
(ion  drifts)  imposed  from  the  magnetosphere,  which  are 
particularly  high  during  Bz  negative  periods. 

We  calculated  the  daily  and  hourly  cross-correlations 

between  Pe  and  Vsw.  Pz.  |P|.  ^sw^zt  VswPyz.  Vsw|P|. 

log(VswlBI).  log(VLlBI).  and  log(VswlB|^). 
where  Byz  is  the  vector  average  of  By  and  Bz,  or  the  square 
root  of  the  sum  of  the  squares  of  their  hourly  averages.  For 
hourly  data,  the  best  correlations  between  Pe  and  most  of 
the  solar  wind  parameters  occurred  when  Pe  lagged  the 
solar  wind  parameters  by  1  h.  Cross-correlations  were 
found  for  Bz  negative,  Bz>=  0,  and  all  Bz  conditions,  for 
each  type  of  solar  wind  structure  for  1978-2008  where  all 
parameters  were  available.  All  cross-correlations  were 
significant  (>0.5)  for  at  least  one  type  of  solar  wind 
structure  for  Bz  negative  or  positive  conditions.  Usually, 
hourly  cross-correlations  were  less  than  daily  cross¬ 
correlations  by  '-0.20  for  Vsw|B|  combinations,  '-0.15  for 
Vsw  or  |B|  alone,  less  than  0.10  for  Bz  combinations,  and 
about  the  same  for  Bz  alone.  The  cross-correlations  were 
larger  for  the  combination  of  the  HSS  and  slow-speed 
solar  wind,  rather  than  for  HSS  or  slow-speed  winds  alone 
for  Vsw  or  combinations  with  Vsw  ^nd  |B|  or  Byz. 

HPe  increases  strongly  with  increasing  Bz  negative 
values  and  with  increasing  Vsw  magnitudes  (Brautigam 
et  al.,  1991;  Liou  et  al.,  1998;  Emery  et  al.,  2008),  so 
correlations  of  Pe  with  Bz  and  Vsw  were  examined  first  as  a 
function  of  solar  wind  structure.  Transients  and  the  slow- 
speed  wind  Pe  correlated  well  with  hourly  Bz  for  Bz 
negative  conditions  with  cross-correlation  coefficients  of 
-0.63  and  -0.64,  respectively,  while  HSS  Pe  cross¬ 
correlations  were  -0.51.  Thus  Bz  negative  is  more  effective 
in  producing  Pe  during  transient  or  slow-speed  wind 
conditions  than  during  HSS.  However,  HSS  and  slow-speed 
wind  Pe  cross-correlations  with  hourly  Vsw  were  0.54 
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compared  to  transient  Pe  cross-correlations  of  0.45, 
showing  that  Vsw  is  relatively  more  effective  for  HSS 
conditions  compared  to  transient  conditions.  The  correla¬ 
tion  of  Pe  with  Vsw  for  positive  Bz  conditions  is  expected 
because  previous  studies  found  HPe  increases  only  slowly 
with  increasing  Bz  positive  while  HPe  doubles  between  Vsw 
of  300  and  600  km/s  (Brautigam  et  al.,  1991 ;  Emery  et  al., 
2008).  Cross-correlations  of  Pe  were  usually  insignificant 
{ <0.5)  for  Bz  positive  conditions  for  any  combination  with 
Bz.  However,  daily  cross-correlations  of  Pe  with  |B|  during 
transient  or  slow-speed  wind  conditions  were  just  over 
0.50,  and  VswBxy  d^ily  cross-correlations  with  Pe  were 
also  significant  for  Bz  positive  and  especially  Bz  negative 
conditions. 

The  best  overall  cross-correlations  of  hourly  Pe  (lagged 
1  h)  and  daily  Pe  for  Bz  positive  and  negative  conditions 
were  found  for  log(V5^|B|),  V^^IBI,  and  Vsw|B|.  The  range  of 
hourly  cross-correlations  of  Pe  with  these  combinations 
are  0.56-0.61, 0.50-0.60,  and  0.56-0.62,  respectively  for  Bz 
negative  conditions  and  0.49-0.57,  0.47-0.60,  and 
0.47-0.53,  respectively  for  Bz  positive  conditions,  where 
the  lower  numbers  are  during  transient  solar  wind 
conditions  and  the  higher  numbers  are  for  the  combina¬ 
tion  of  HSS  and  slow-speed  wind.  V^^B^  is  expected  from 
theoretical  considerations  to  be  a  better  choice  for  a 
power  input  like  Pe  (e.g.  Table  2  of  Gonzalez  et  al.,  1994) 
because  it  it  is  part  of  the  merging  energy  transfer  term  at 
the  magnetopause.  However,  the  interplanetary  electric 
field  Vsw|B|  also  correlates  well  with  Pe  as  with  y^p.  |B| 
combinations  in  particular  correlate  well  for  HSS  because 
the  geomagnetic  disturbances  are  caused  by  the  south¬ 
ward  component  of  large-amplitude  Alfven  waves  in  the 
solar  wind  (Tsurutani  et  al.,  1995),  and  so  the  fluctuations 
in  the  IMF  and  thus  |B|  are  more  important  for  HSS 
compared  to  the  relatively  slow  variations  in  IMF 
experienced  in  transient  solar  structures. 

Fig.  3a  and  b  show  the  daily  average  Vsw|B|  values  in 
mV/m  against  the  daily  average  Pe  values  between  1978 
and  2008  for  Bz  negative  and  positive  conditions.  The  daily 
(solid  line)  and  hourly  (dashed  line)  fits  are  also  plotted, 
where  the  HSS  and  slow-speed  wind  fits  are  linear  and  the 
transient  fits  are  quadratic.  The  Pe  values  are  larger  for  Bz 
negative  conditions,  as  expected,  where  the  largest  Vsw|B| 
values  are  during  transients  with  large  |B|  values  as  shown 
in  Fig.  2b.  Usually  the  largest  Pe  values  are  brought  by 
transient  solar  wind  structures,  although  HSS  in  Bz 
positive  conditions  can  precipitate  relatively  large  Pe. 
The  hourly  fits  are  for  larger  values  of  Pe  than  the  daily  fits 
for  Bz  negative  conditions,  while  the  hourly  fits  are  for 
smaller  values  of  Pe  for  Bz  positive  conditions.  This  shows 
that  an  equivalent  increase  in  the  IMF  magnitude  (or  B^) 
results  in  larger  Pe  for  Bz  negative  conditions  and  not 
much  increase  in  Pe  for  Bz  positive  conditions,  as  expected. 

For  Bz  positive  conditions,  the  lower  Pe  values  are 
during  transient  solar  wind  structures,  showing  that  |B| 
increases  are  not  as  effective  for  Pe  during  transient 
structures  compared  to  HSS  structures.  For  Bz  negative 
conditions,  the  difference  between  the  Pe  during  transient 
structures  is  a  little  higher  compared  to  slow-speed  and 
HSS  structures  for  interplanetary  electric  fields  less  than 
6mV/m.  indicating  that  the  same  |B|  (or  Bz)  is  more 
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Fig.  3.  (a)  The  daily  averages  of  the  hourly  products  of  solar  wind  speed 
(Vsw)  and  scalar  IMF  |5|  in  three  types  of  daily  solar  wind  structures 
versus  the  daily  average  auroral  electron  power  (Pe)  for  average  daily 
IMF  negative  conditions  from  November  1978  to  August  2008  for  all 
hours  with  Vjw,  IMF,  and  Pg.  Pe  lags  Vsw  and  IMF  by  1  h.  The  green  lines 
are  linear  fits  for  HSS  combined  with  slow-speed  winds  for  the  plotted 
daily  averages  (solid)  and  for  the  hourly  fit  (dashed).  The  orange  lines 
are  quadratic  fits  of  Pe  to  Vsw|6|  for  daily  (solid)  and  hourly  (dashed) 
averages  with  transient  solar  structures.  The  corresponding  daily  and 
hourly  cross-correlations  are  also  listed,  (b)  Same  as  (a)  for  Bz  positive 
conditions. 


effective  in  depositing  larger  Pe  for  transients,  probably 
because  the  Bz  negative  conditions  last  longer.  The 
difference  in  transient  Pe  magnitudes  compared  to 
HSS  Pe  magnitudes  as  a  function  of  the  sign  of  Bz 
emphasizes  the  different  characteristics  of  transients 
compared  to  HSS. 

For  Bz  negative  conditions,  slightly  better  cross-corre¬ 
lations  with  hourly  Pe  were  found  with  V^^Bz  (r  =  -0.64) 
and  VswBz  (r=-0.62)  for  HSS  and  slow-speed  wind 
conditions,  and  much  better  cross-correlations  were 
found  with  VswBz  {r=-0.65).  B,  {r=-0.63).  Vl^Bz 
(r  =  -0.60),  and  VswByz  (r  =  0.58)  during  hours  with 
transient  solar  structures.  These  are  similar  to  the  hourly 
cross-correlations  of  Vs^Bz  and  V^^Bz  with  HPe  in  the 
southern  hemisphere  found  in  Emery  et  al.  (2008). 


Daily  Averages  1 978-2008 


transierit 

OHSSard  O  jlov  _daiiy(r«0  80X --t>ourV(r.0  62)  ^ 

6  i  6 


5  10 

VswlBI  (mV/m)  for  B2<0 


5  10 

Vsw|B|  (mV/m)  for  B2>-0 


BA.  Emery  et  at.  /  joumol  of  Atmospheric  ond  Solar-Terrestrial  Physics  71  (2009)  1157-1175 


1167 


Because  Pc  is  larger  for  negative  conditions  than  for  Bz 
positive  conditions,  the  cross-correlations  of  hourly  Pe 
with  all  Bz  conditions  is  often  worse  than  the  cross¬ 
correlations  for  either  Bz  condition  because  of  the 
increased  scatter  in  the  Pe  values.  Bz  positive  cross¬ 
correlations  for  Pc  are  largest  for  HSS  structures  compared 
to  transient  structures,  while  Bz  negative  cross-correla¬ 
tions  are  similar  for  combinations  including  Bz. 

Charan  Dwivedi  et  al.  (2009)  found  cross-correlations 
of  daily  Ap  with  V4w|P|  to  be  relatively  small  in  even  SC 
ascending  phases  (~0.70  in  A20  and  '--0.65  in  A22)  and 
relatively  large  in  odd  SC  ascending  phases  (~0.75  in  A21 
and  ~0.80  in  A23).  We  found  the  cross-correlations  of 
daily  Pc  with  Vsw|P|  were  0.56  in  1987-1988  (~A22),  and 
were  0.73  in  1997-1998  (~A23),  for  all  Bz  and  solar  wind 
structure  conditions.  The  transient  cross-correlations  of  Pe 
with  Vsw|B|  were  0.07  worse  than  those  for  HSS  and  slow- 
speed  winds  in  A22,  and  better  than  them  by  0.05  in  A23 
when  the  percentage  of  transients  nearly  doubled.  How¬ 
ever,  part  of  the  reason  for  the  poor  correlations  in 
1987-1988  could  be  that  the  solar  wind  plasma  or  IMF 
observations  were  missing  on  33%  of  the  days. 

4.  Periodic  analyses  of  solar  wind  and  auroral 
electron  power 

4.1.  Lomb-Scargle  analyses 

Several  authors  have  analyzed  the  spectral  character¬ 
istics  of  the  solar  wind.  Svalgaard  and  Wilcox  (1975) 
studied  the  IMF  Bx  periodicity  over  five  solar  cycles  and 
found  a  strong  27-day  periodicity.  Gosling  et  al.  (1976, 
1977)  reported  27-day  periodicities  in  Vsw,  with  peaks  at 
13.5  days  for  certain  years,  and  weak  or  absent  7-day 
periodicities.  Similarly,  several  authors  found  27-,  13.5- 
day  (one  half  of  27-day),  9-day  (one  third  of  27-day)  and 
7-day  (one  fourth  of  27-day)  periodicities  in  the  IMF  and 
Vsw  (Neugebauer  et  al.,  2000;  Nayar  et  al.,  2001;  Verma, 
2001;  Temmer  et  al.,  2007).  Mursula  and  Zieger  (1996) 
showed  the  27-  and  13.5-day  periodicities  in  Vsw.  IMF 
sector  component  and  Kp,  and  ascribed  the  13.5-day 
periodicity  to  the  occurrence  at  1  AU  of  two  HSS  per  solar 
rotation  corresponding  to  two  different  magnetic  polarity 
sectors  with  a  tilted  heliosheet.  Lindblad  (1981)  showed  a 
preference  for  HSS  on  Bartels  day  17  for  Bx  negative 
(toward  the  Sun),  and  a  preference  for  Bartels  day  4  for  Bx 
positive  (away  from  the  Sun),  corresponding  to  two  HSS 
per  solar  rotation.  Rangarajan  and  Mavromichalaki  (1989) 
found  similar  preferences  for  HSS.  with  additional  pre¬ 
ferences  for  transients  on  Bartels  days  11  and  25  for  Bx 
negative  and  positive  sectors. 

Because  of  the  high  correlations  of  Pe  and  HPe  with  3-h 
Kp  and  Op  (Section  3.4  and  Emery  et  al.,  2008),  we 
anticipate  similar  periodicities  in  geomagnetic  indices  and 
Pe.  The  major  periodicities  of  the  Earth’s  geomagnetic 
indices  Kp,  Ap,  or  aa  are  for  11.1 -year,  17.5,  12.7,  and  6 
months,  with  shorter  27-  and  13.5-day  periods  (e.g. 
Delouis  and  Mayaud,  1975;  Gonzalez  et  al.,  1990;  Gorney, 
1990;  Prestes  et  al.,  2006).  The  22-year  Hale  cycle  is  seen 
in  some  studies,  and  not  seen  in  others  (Prestes  et  al.. 


2006).  The  6-month  periodicity  is  quite  marked,  especially 
after  1930  when  the  solar  variability  increased 
during  solar  minima  (Prestes  et  al.,  2006).  Oliver  et  al. 
(2000)  attributed  ^^65%  of  this  semiannual  variation 
to  the  equinoctial  effect  based  on  the  angle  between  the 
solar  wind  flow  direction  and  the  Earth’s  dipole  axis  that 
makes  the  IMF  coupling  ~25%  less  effective  at  the 
solstices.  The  remaining  6-month  periodicity  was  attrib¬ 
uted  ~20%  to  the  axial  effect  and  --15%  to  the  Russell- 
McPherron  effect  (Russell  and  McPherron,  1973).  The  axial 
effect  reflects  changes  in  the  solar  wind  parameters 
measured  at  1  AU  because  the  Earth  is  at  the  highest 
heliospheric  latitudes  on  6  March  (-7.2  )  and  on  6 
September  (+7.2’)  (e.g.  Oliver  et  al.,  1996),  where  it  is 
more  in  line  with  solar  midlatitude  coronal  holes  (e.g. 
Bohlin,  1977).  Oliver  et  al.  (2000)  found  a  6km/s 
semiannual  variation  in  Vsw  from  1963  to  1997  on  a 
background  average  Vsw  of  — 440km/s. 

From  the  previous  sections,  we  anticipate  that  some  of 
the  periodicities  in  Vjw  will  be  in  Pe.  We  also  anticipate 
finding  an  enhanced  semiannual  periodicity  (—183  days) 
similar  to  that  found  in  the  Kp  index  for  Pg,  which 
maximizes  in  the  equinoxes  when  magnetospheric  re¬ 
connection  is  increased  compared  with  solstices  (e.g. 
Cliver  et  al.,  2000). 

We  used  a  Lomb-Scargle  analysis  (Lomb,  1976;  Scargle, 
1982)  of  hourly  total  Vsw  and  Pp  because  a  Lomb  analysis 
can  handle  the  significant  gaps  in  Vsw  data  before  about 
1995  (the  launch  of  the  WIND  spacecraft).  The  full  Pe  data 
were  used,  since  there  were  very  few  data  gaps.  Several 
authors  have  used  and  discussed  the  advantages  of  the 
Lomb  analysis  over  the  traditional  fast  Fourier  transform 
(FFT)  method  or  other  methods  (e.g.  Harrington  et  al., 
1996;  Vityazev,  1997;  Laguna  et  al.,  1998;  Thong  et  al.. 
2004).  Thong  et  al.  (2004)  averaged  different  Lomb- 
Scargle  transforms  after  de-normalization  of  the  Lomb 
transform,  which  is  variance  (cr^)  normalized  in  spectral 
power  y.  The  spectral  power  y  is  converted  to  amplitude  A 
as  A  =  2y/(ya^/N),  where  N  is  the  number  of  sample 
points. 

We  want  to  compare  the  hourly  Vsw  3nd  Pe  amplitudes 
at  several  periodicities  over  several  solar  cycles.  To  test 
our  Lomb-Scargle  analyses,  we  calculated  FFT  amplitudes 
using  daily  Vsw  ^nd  Pe  values  for  each  year  of  data  where 
data  gaps  were  filled  in.  The  amplitudes  and  periods  from 
both  Lomb-Scargle  and  FFT  analyses  were  similar, 
although  the  Lomb  analyses  had  better  frequency  resolu¬ 
tion.  The  Lomb  results  were  also  similar  using  hourly 
values  in  place  of  daily  values.  Thus,  the  following 
Lomb-Scargle  periodic  analyses  use  hourly  Vsw  ^nd  Pe 
over  various  time  intervals. 

We  are  interested  in  the  periodicities  and  in  the 
maximum  amplitudes  of  these  periodicities.  Figs  4a-d 
show  the  Lomb  analysis  of  the  hourly  Vsw  and  Pe  in  the 
descending  phase  of  SC  23  from  02190  to  06027  when  HSS 
were  present  for  50%  of  the  time,  with  transients  present 
16%  of  the  time  and  slow-speed  wind  the  remaining  34% 
of  the  time,  as  shown  in  Fig.  la.  Fig.  4a  shows  the  total  Vsw 
and  Pe,  while  Fig.  4b  is  only  for  the  hours  with  transients. 
Fig.  4c  for  the  HSS  hours,  and  Fig.  4d  for  the  slow-speed 
hours.  As  anticipated,  the  transient  periodicities  in  Fig.  4b 
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Period  (days)  from  hourly  data  Period  (days)  from  hourly  data 


Fig.  4.  (a)  Lomb-Scargle  spectral  amplitudes  using  all  the  hourly  V^sw  and  P*  values  during  SC  23  descending  phase  (02190-06027).  (b)  Same  as  (a)  for  all 
haurs  with  transients,  (c)  Same  as  (a)  for  all  hours  with  HSS.  (d)  Same  as  (a)  far  all  hours  with  slaw  speed  wind. 


are  much  noisier  than  any  of  the  other  solar  structures. 
The  periodicities  around  5,  7,  9, 13.5,  and  27  days  seen  in 
all  Vsw  3nd  Pe  data  in  Fig.  4a  are  clearest  in  Fig.  4c  for  HSS 
hours,  suggesting  that  HSS  structures  impose  the  most 
structure  on  Vsw  ^nd  Pp.  These  periodicities  are  also  seen 
in  |B|,  Bx,  and  GSM  By,  but  are  not  evident  in  GSM  B^  (not 
shown).  For  B^  and  there  are  strong  amplitudes  of  13.5 
and  27  days  reflecting  the  sector  structure  of  a  tilted 
heliosheet  (e.g.  Mursula  and  Zieger,  1996).  The  HSS  and 
total  spectra  are  especially  close  for  Vsw.  while  for  Pe,  7-, 
9-  and  25-day  peaks  are  also  relatively  large  in  the  slow- 
speed  wind.  The  Pe  peak  at  --208  days  in  Fig.  4a  appears 
to  be  shifted  from  ‘semiannual’  (~183  days)  peaks  in  the 
HSS  and  slow-speed  hours  in  Figs.  4c  and  d  by  the 
significant  amplitudes  ^-200  days  seen  in  the  transient 
hours  in  Fig.  4b. 


We  made  similar  analyses  for  other  solar  cycle  phase 
intervals  indicated  in  Table  1.  In  general,  the  7-,  9-,  13.5-, 
25-,  27-,  and  29-day  periodicities  in  Vjw  usually  also  appear 
in  Pp.  The  longer  periodicities  are  sometimes  in  both,  and 
sometimes  are  not.  Figs.  5a  and  b  are  the  periodicities  for 
the  total  Vsw  ^nd  Pp  for  366  days  in  the  SC  22/23  minimum 
from  96031  to  97030,  and  in  the  SC  23/24  minimum  from 
07209  to  08209,  respectively.  Gibson  et  al.  (2009)  discuss 
these  solar  minimum  periodicities  in  the  context  of  the 
solar  drivers  and  their  effects  on  the  whole  geospace 
system.  Both  solar  minima  show  periodicities  in  the  total 
Vsw  3nd  Pp  around  27,  13.5,  9,  7  and  5  days,  but  the 
amplitudes  in  SC  23/24  are  much  stronger  and  cleaner.  In  SC 
minimum  22/23,  there  is  also  a  very  strong  ‘semiannual’ 
peak  in  both  Vsw  ^nd  Pp.  We  will  discuss  the  differences  in 
these  strong  solar  minimum  amplitudes  in  Section  4.3. 


BA  Emery  ec  oi  /  Joumol  of  Atmospheric  ond  Solar-Terrestriol  Physics  71  (2009)  1157- 1175 


1169 


3  Solar  Minimum  96031  -97030 


Period  (days)  from  hourly  data 


b  Solar  Minimum  07209-08209 


Fig.  5.  (a)  Lomb-Scargle  spectral  amplitudes  using  all  the  hourly  VTs^^and 
Pc  values  in  one  year  during  SC  22/23  minimum  (96031-97030).  (b) 
Same  as  (a)  for  one  year  during  SC  23/24  minimum  (07209-08209). 


4.2.  Solar  rotation  penodidf/es 

Fig.  6a  shows  the  Lomb  amplitudes  for  the  total  hourly 
Vsw  red  triangles  and  Pe  as  dark  green  x’s  in  the  solar 
rotation  and  related  periodicities  of  7-,  9-,  13.5-,  27-day, 
and  the  maximum  of  the  25-  and  29-day  amplitudes. 
These  periodicities  are  almost  always  present  in  \4w  and 
in  Pe,  and  are  correlated.  The  Lomb  amplitudes  using 
hourly  Vsw  and  Pe  values  over  yearly  intervals  are  plotted 
as  magenta  diamonds  and  cyan  circles.  The  y-axes  are  in 
km/s  (Vsw)  and  GW  (Pe),  where  a  dotted  line  is  drawn  at 
5  GW  (~33.8  km/s)  as  guidance. 

Cliver  et  al.  (1996)  studied  the  daily  and  longer  aa 
index  using  fast  Fourier  transform  (FFT)  analyses  for  six 
Hale  cycles  from  1844  to  1994.  They  found  the  aa  activity 
to  be  enhanced  in  the  descending  and  minimum  (DN) 
even  cycles  and  in  the  ascending  and  maximum  (AX)  odd 
cycles  in  the  middle  of  each  Hale  cycle.  They  showed  that 


the  182-,  27-,  and  13.5-day  amplitudes  were  largest  for  DN 
even  cycles,  with  9-day  amplitudes  about  the  same  for  DN 
even  and  odd  cycles.  They  suggested,  that  in  addition  to 
the  enhanced  GSM  Bz  negative  from  the  Russell-McPher- 
ron  and  Rosenberg-Coleman  effects,  that  faster  and  longer 
duration  HSS  in  DN  even  cycles  accounted  for  the  DN 
increases  in  aa,  and  that  CMEs  in  AX  odd  cycles  produced 
the  aa  increases  in  odd  AX  cycles.  Cliver  et  al.  (1996)  also 
found  that  the  adjacent  25-  and  29-day  periodicities  in  aa 
were  larger  and  more  common  for  DN  even  cycles.  Again, 
because  the  22-23  Hale  cycle  is  unusual,  these  findings 
may  not  hold  for  in  this  Hale  cycle. 

The  ‘importance’  parameter  of  Maris  and  Maris  (2005) 
was  highest  for  DN20  followed  by  DN21  for  HSS,  which 
confirms  the  Cliver  et  al.  (1996)  results  of  more  active  HSS 
in  DN  even  cycles.  But  for  transients,  the  ‘importance’ 
parameter  in  Maris  and  Maris  (2005)  was  highest  in  AX20, 
followed  by  AX22  and  finally  by  AX21,  which  is  against  the 
Cliver  et  al.  (1996)  speculation  that  transients  contribute 
most  in  AX  phases  of  odd  cycles. 

From  our  study,  we  see  in  Fig.  6a  that  the  27-day  and 
the  maximum  of  the  subsidiary  25-  and  29-day  periodi¬ 
cities  are  relatively  strong  for  Vsw  and  Pe  in  all  DN  cycles, 
especially  in  the  yearly  amplitudes,  which  is  different 
from  the  aa  index  findings  of  Cliver  et  al.  (1996).  Similarly, 
the  7-,  9-  and  13.5-day  periodicities  are  strongest  in  all  DN 
phases.  The  strongest  13.5-day  periodicity  in  Vjw  is  during 
DN20  in  agreement  with  Cliver  et  al.  (1996). 

We  also  see  in  Fig.  6a  strong  9-day  periodicities  in  Vsw 
and  Pe  during  DN23,  especially  in  2005.  This  is  in 
agreement  with  the  observations  of  Temmer  et  al. 
(2007),  who  found  the  9-day  period  in  2005  was  caused 
by  three  CHs  distributed  ~120^  apart  in  solar  longitude. 
Mlynczak  et  al.  (2008)  found  a  similar  9-day  periodicity  in 
the  infrared  cooling  from  2002  to  2006  by  [NO]  and  [CO2I 
between  100  and  200  km.  They  found  a  9-day  periodicity 
in  daily  Ap  and  iCp,  but  did  not  find  any  9-day  periodicity  in 
the  solar  ultraviolet  (UV)  data.  Lei  et  al.  (2008a.b)  and 
Thayer  et  al.  (2008)  also  found  9-  and  7-day  periodicities 
in  the  neutral  density  ~400  km  from  the  CHAMP  satellite 
between  2004  and  2007.  The  7-  and  9-day  periodicities  in 
Vsw  3nd  Pe  were  also  strong  for  DN23  as  shown  in  Figs.  4a, 
5b.  and  6a.  In  Figs.  2a  and  c,  we  see  exceptionally  large 
contributions  from  HSS  to  the  total  and  Pe  in  N23/24. 
There  are  correspondingly  lower  contributions  to  the  total 
Vsw  and  Pe  by  the  slow  solar  wind. 


4.3.  Semiannual  and  other  longer  periodicities 

In  Fig.  6b,  we  plot  Lomb  amplitudes  in  km/s  and  GW  of 
periodicities  longer  than  100  days  in  the  hourly  Vsw  and  Pe 
during  SC  phases  and  yearly  intervals.  The  yearly 
amplitudes  are  shown  in  the  top  panel.  The  longest 
periodicity  shown  is  for  the  1.3-1.4  year  (465-520  days  or 
15.2-17.1  months)  periodicity  in  the  Sun  that  was 
observed  in  the  solar  wind  starting  in  1987  by  Richardson 
et  al.  (1994).  In  Fig.  6b,  we  see  it  relatively  strong  in  Vsw 
and  Pe  for  SC  22,  but  not  always  in  the  same  phases. 
Prestes  et  al.  (2006)  have  a  relatively  strong  periodicity  in 
monthly  aa  and  Ap  at  15.9  and  16.9  months  in  their  Fig.  5. 
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3  Amplitude  of  Periodicities 


b  Amplitude  of  Periodicities 


Fig-  6.  (a)  Amplitudes  from  the  Lomb-Scargle  analyses  of  all  hourly  P*  (as  green  x's)  and  all  hourly  (as  red  triangles)  for  7-.  9-,  13.5-,  27-day.  and  the 
maximum  of  25  or  29-day  periods  related  to  the  solar  rotation  period  over  the  descending  (D),  minimum  (N),  ascending  (A)  and  maximum  (X)  phases  of 
the  solar  cycle.  Yearly  Lomb-Scargle  amplitudes  of  Vsw  are  magenta  diamonds  and  are  cyan  circles  for  P^  at  the  midpoint  of  each  year,  (b)  Same  as  (a)  for 
periods  of  120-145  days,  166-200  days  (semiannual),  270-300  days,  and  465-520  days  (1.3-1.4  years)  from  the  hourly  Lomb-Scargle  analyses.  Yearly 
Lomb-Scargle  amplitudes  as  diamonds  (Vsw)  and  circles  (P^)  are  dashed  lines  in  the  top  panel.  Yearly  Lomb-Scargle  amplitudes  between  145  and  165  days 
are  shown  as  inverted  triangles  ( Vsw)  and  squares  (Pe)  in  the  plot  for  120-145  days.  These  same  symbols  are  used  with  dashed  lines  in  the  second  panel  for 
equinox  semiannual  (sa)  amplitudes  from  polynomial  fits  to  daily  Vsw  and  Pe  data  sorted  by  daynumber.  Solstice  semiannual  amplitudes  in  are  marked 
with  ‘S'. 


Average  Daily  Auroral  Electron  Power  (GW)  CT  Average  Daily  Solar  Wind  Velocity  (km/s) 
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Fig.  7.  (a)  Average  daily  solar  wind  velocity  from  95109  to  99270  (SC  22/23  minimum  and  SC  23  ascending  phase)  showing  the  averages  of  the  FFT  and 
Lomb-Scargle  semiannual  fits  in  blue,  and  the  27-day  fits  in  magenta  on  top  of  annual  averages,  with  several  CMEs  labeled,  (b)  Same  as  (a)  for  the  daily 
average  auroral  electron  power. 
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Silverman  and  Shapiro  (1983)  showed  this  1.4-year 
periodicity  in  auroral  sightings  in  Sweden  from  1721  to 
1943,  with  peaks  approximately  every  65-68  years 
(~1764,  M832,  ~1878,  and  by  extrapolation  ~1948  and 
~2013).  The  reason  for  this  1.3-1.4-year  variability  is  not 
known. 

The  periodicity  for  270-300  days  (8.9-9.9  months) 
plotted  in  the  third  panel  of  Fig.  6b  is  relatively  common 
in  the  hourly  Vsw.  with  significant  amplitudes  at  times. 
This  periodicity  is  only  weakly  present  in  Pe-  Prestes  et  al. 
(2006)  show  a  ~9.1 -month  periodicity  in  monthly  aa  and 
Ap  in  their  Fig.  5. 

Richardson  and  Cane  (2005)  described  an  intermittent 
~1 50-day  quasi-periodicity  in  interplanetary  and  solar 
phenomena  during  SC  23.  Lean  (1990)  noted  that  it 
usually  appeared  around  the  maxima  of  cycles  12-21  in 
sunspot  areas  with  a  period  between  130  and  185  days. 
We  split  this  ~1 50-day  periodicity  into  three  sections 
between  120  and  145  days  (3.9-4.8  months)  and  the 
semiannual  periodicity  between  166  and  200  days 
(5. 5-6.6  months),  with  the  intervening  periodicity  of 
145-165  days  (4.8-5.5  months)  plotted  as  inverted 
triangles  (Vjw)  and  squares  (Pe)  on  the  145-165  day  plot 
in  Fig.  6b  in  the  top  panel.  Periodicities  in  V^sw  and  Pe 
between  120-145  days  and  145-165  days  are  relatively 
weak  compared  with  the  solar  rotation  periodicities  in 
Fig.  6a.  These  periodicities  were  listed  by  Prestes  et  al. 
(2006)  for  monthly  aa  and  Ap  from  1932  to  2001  as  4.7 
and  4.9  months  (143  and  149  days). 

The  ’semiannual’  amplitudes  are  shown  in  the  second 
panel  of  Fig.  6b.  The  largest  ’semiannual’  Lomb  power  in 
Vsw  is  for  N21/22  at  56km/s  for  189  days.  The  largest  Pe 
amplitudes  are  yearly  peaks  in  1982  at  174  days,  in  1996  at 
193  days,  and  in  1999  at  192  days.  The  ’semiannual’  peaks 
in  Vsw  and  Pe  during  96031-97030  are  also  shown  in 
Fig.  5a.  Compared  with  amplitudes  at  other  periodicities, 
the  yearly  ‘semiannual’  amplitudes  in  Pe  are  relatively 
large  compared  with  the  Usw  amplitudes  in  Figs.  5a 
and  6b. 

Because  the  Lomb-Scargle  analysis  does  not  specify 
the  phase  of  the  amplitude  peaks,  semiannual  polynomial 
fits  were  made  for  all  daily  Vsw  and  Pe  data  sorted  in 
daynumber.  For  Vjw  between  1972  and  2008,  we  found  a 
semiannual  amplitude  of  5.7km/s  for  equinoctial 
peaks  (22  March  and  22  September)  and  an  amplitude 
of  5.1  km/s  for  solar  axial  peaks  (6  March  and  6 
September)  with  an  average  Vsw  of  447  km/s  or  ~2.6% 
equinoctial  variation.  These  values  are  consistent  with  the 
6  km/s  equinoctial  amplitude  found  by  Cliver  et  al.  (2000) 
for  V^sw  between  1963  and  1997.  The  semiannaul  equinoc¬ 
tial  amplitude  for  all  Pe  data  was  2.7  GW  with  an  average 
Pe  of  34.5  GW  or  --15.7%  variation. 

Semiannual  polynomial  fits  to  Vsw  and  Pe  were  made 
for  all  SC  phases.  The  equinoctial  semiannual  amplitudes 
are  plotted  in  the  second  panel  of  Fig.  6b,  along 
with  solstice  semiannual  amplitudes  that  were  found  for 
D22,  D23  and  N23/24.  The  average  Vsw  equinoctial 
amplitude  from  the  seven  SC  phases  with  equinox  peaks 
over  the  entire  37  years  was  ~6.8  km/s  (—3.0%  variation), 
similar  to  the  5.7  km/s  found  using  all  Vsw  d^ta.  Pe 
had  equinoctial  polynomial  peaks  for  all  SC  phases 


except  D23,  where  the  solstice  Vsw  variations  cancelled 
out  the  normal  equinoctial  peaks  in  Pe  due  to  preferred 
magnetospheric  reconnection  (Cliver  et  al.,  2000).  The 
average  equinoctial  amplitude  in  Pe  from  the  other  SC 
phases  was  —2.6  GW,  close  to  the  2.7  GW  found  using 
all  Pe  data. 

Fig.  7a  and  b  show  the  daily  ^nd  Pe  during  the 
NA23  period  from  95109  to  99270.  The  Lomb  and  FFT 
estimated  semiannual  and  27-day  amplitudes  were  very 
similar  (as  expected)  so  they  were  averaged  together  and 
plotted  on  top  of  annual  averages.  The  lowest  annual  Vsw 
values  were  in  1997,  while  the  lowest  Pe  values  were  in 
1996  and  1997.  The  27-day  variation  was  set  to  agree  with 
the  Vsw  27-day  variation  at  the  start  of  the  minimum 
period,  and  so  is  not  correct  in  the  ascending  phase, 
especially  when  strong  CMEs  were  present,  or  when  a 
strong  13.5-day  periodicity  shifted  the  fit  to  the  27-day 
periodicity  half  a  cycle. 

The  semiannual  variation  was  timed  to  maximize  in 
the  equinoxes  since  we  calculated  polynomial  equinoctial 
peaks  for  both  Vsw  ('^IQ  km/s)  and  Pe  (—5.0 GW)  in  NA23, 
which  are  similar  to  the  Lomb-Scargle  semiannual 
amplitudes  in  Fig.  6b  of  —20  km/s  and  —5.3  GW  used  in 
Figs.  7a  and  b.  Watari  and  Watanabe  (2000)  remarked  on 
the  low  geomagnetic  activity  during  solstices  between 
1995  and  1997,  with  very  low  values  in  the  solstice 
of  1996.5,  which  is  also  evident  in  Pe  and  Vsw  in  Figs.  7b 
and  a. 


5.  Conclusions 

We  have  shown  that  the  solar  wind  structure  plays  a 
determining  role  in  the  magnitude  of  the  electron  auroral 
input  as  characterized  by  the  hourly  average  auroral 
electron  power  (Pe)  over  three  solar  cycles  (SCs).  The 
variability  in  the  total  magnitude  of  Pe  follows  the 
variations  in  Pe  during  periods  of  high-speed  streams 
(HSS).  This  is  consistent  with  the  similar  findings  by 
Richardson  et  al.  (2002)  that  the  variations  in  the  averages 
of  the  3-h  aa  index,  the  IMF  magnitude,  and  the  solar  wind 
velocity  closely  follow  the  variations  associated  only  with 
HSS.  Many  of  the  Vsw  solar  periodicities  in  HSS  are 
transferred  to  Pe  more  effectively  than  in  slow-speed  solar 
wind  periods,  or  in  times  of  transients  when  the 
periodicities  are  noisy.  This  finding  about  the  dominance 
of  HSS  in  time  variations  is  not  too  surprising  since  HSS 
occur  43%  of  the  time,  and  the  highest  velocities,  by 
definition,  are  during  HSS.  Pe  and  aa  both  increase  with 
increasing  Vsw*  although  correlations  with  negative  are 
generally  more  important. 

The  characteristics  of  the  hourly  average  solar  wind 
velocity  (Vsw)  ^nd  Pe  are  similar  over  the  SCs.  HSS  from 
coronal  holes  (CH)  contribute  —47%  to  the  total  Pe  and  Vsw. 
and  delineate  most  of  the  structures  seen  in  both  with 
large  27-day  cross-correlations  of  0.8  for  both.  The  largest 
HSS  contributions  to  Pe  are  —57%  in  the  descending  phases 
of  the  SC,  with  the  smallest  contributions  of  —32%  in  solar 
maxima.  The  slow-speed  wind  structures  contribute  —33% 
to  Pe  in  solar  minima,  and  —24%  in  solar  maximum  and 
descending  phases. 
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The  contribution  of  transient  (CME)  solar  wind  struc¬ 
tures  is  much  more  variable.  Transients  make  their 
maximum  contribution  during  solar  maxima  (X),  con¬ 
tributing  ~43%  to  Pe  in  X21  and  X22,  and  ^38%  in  X23 
maximum.  During  solar  minimums,  transients  only  con¬ 
tribute  ~6%  to  Pe. 

The  cross-correlations  of  Vsw|B|  in  slow-speed  wind 
structures  together  with  HSS  with  hourly  average  Pe 
lagging  by  1  h  are  0.53  and  0.62  for  positive  and 
negative  conditions,  respectively,  and  can  be  represented 
by  linear  equations.  Because  transients  represent  the 
largest  |B|  values,  quadratic  fits  are  better,  and  the  cross¬ 
correlations  of  Pe  with  Vswl^l  3re  lower  at  0.47  and  0.56  for 
Bz  positive  and  negative  conditions,  respectively.  Transi¬ 
ents  are  less  effective  than  other  solar  wind  structures  at 
producing  Pe  for  Bz  positive  conditions.  Bz  negative 
conditions  are  relatively  more  important  for  transients, 
and  the  highest  hourly  cross-correlations  of  transient  Pe 
are  -0.65  with  negative  VswBz,  The  cross-correlations  of 
hourly  Pe  with  3-h  Up  and  Kp  are  higher  at  0.68  and  0.72, 
respectively,  for  all  Pe  data  from  1978  to  2008. 

The  solar  rotational  27-,  13.5-,  9-,  and  7-day  periodi¬ 
cities  are  almost  always  present  in  the  hourly  total  Vsw 
and  Pe.  The  magnitudes  are  usually  largest  in  the 
descending  (D)  and  minimum  (N)  phases  of  each  solar 
cycle.  DN23  showed  large  7-  and  9-day  periodicities  in  Vsw 
and  Pe.  These  same  periodicities  also  appeared  in  the 
neutral  density  at  400  km  (Lei  et  al.,  2008a, b;  Thayer  et  al., 
2008),  and  in  the  power  radiated  through  infrared  cooling 
by  ICO2]  and  [NO]  between  100  and  200km  (Mlynczak 
et  al..  2008). 

Relatively  weak  amplitude  120-145  or  145-165  days 
periodicities  are  in  the  total  hourly  and  Pe.  These  could 
be  related  to  the  quasi-periodic  150-day  periodicities  seen 
in  sunspots  and  other  phenomena.  Longer  periodicities  of 
270-300  days  were  relatively  common  for  Vsw.  and  less 
common  for  Pe.  The  1.3- 1.4-year  periodicity  was  found  in 
X22  around  1990  in  Pe  and  Vsw 

The  ‘semiannual’  (183  +  17  day)  amplitudes  in  total 
hourly  Vsw  were  checked  with  equinoctial  polynomial  fits. 
The  Vsw  magnitude  at  1  AU  is  affected  by  the  axial  tilt  of 
the  Earth’s  orbit,  with  the  highest  heliospheric  latitudes 
seen  on  6  March  and  6  September.  Equinoctial  (22  March 
and  22  September)  and  axial  amplitudes  were  similar  for 
Vsw  The  largest  equinoctial  amplitudes  in  Vsw  were  for 
N21/22  (-31  km/s),  N22/23  (-21  km/s)  and  A23  (-16  km/ 
s),  with  correspondingly  high  equinoctial  amplitudes  in  Pe 
of  4.4,  5.1  and  4.9  GW,  respectively.  Four  other  SC  phases 
showed  Vsw  equinoctial  amplitudes.  D22  showed  a  strong 
Vsw  solstice  amplitude  (—25  km/s),  with  weaker  Vsw 
solstice  amplitudes  in  D23  and  N23/24.  Equinoctial  fits 
were  found  in  Pe  for  all  SC  phases  except  for  D23.  The 
average  equinoctial  semiannual  amplitudes  were  —6  km/s 
(—2.6%  variation)  in  Vsw  and  —2.7 GW  (—15.7%  variation) 
in  Pe,  where  the  preferred  magnetospheric  reconnection  at 
equinoxes  (Cliver  et  al.,  2000)  greatly  enhanced  the 
equinoctial  amplitudes  in  Pe. 

The  previous  six  Hale  cycles  from  10-11  to  20-21 
showed  higher  aa  indices  in  DN  even  cycles  and  AX  odd 
cycles,  except  for  Hale  cycle  16-17  from  —1925  to  —1945 
(Cliver  et  al.,  1996).  This  aa  pattern  was  broken  again  in 


the  last  22-23  Hale  cycle  (Kozyra  et  al.,  2006)  and  was 
relatively  weak  in  Hale  cycle  20-21  (Fig.  4  of  Cliver  et  al., 
1996).  We  found  the  total  Vsw  3nd  Pe  values  to  be  usually 
largest  in  all  D  phases  because  of  the  dominance  of  HSS  at 
this  time.  There  was  an  increase  of  HSS  in  N23/24 
compared  with  N22/23  in  the  abnormal  Hale  cycle,  where 
the  total  and  HSS  Vsw  magnitudes  were  increased  —9% 
from  N22/23.  However,  the  Pe  magnitudes  were  about  the 
same  because  the  enhanced  Vsw  was  offset  by  a  decrease 
in  IMF  magnitude  of —14%  in  N23/24  compared  to  N22/23. 

Cliver  et  al.  (1996)  found  the  semiannual,  13.5-,  27-day, 
and  the  subsidiary  25-  and  29-day  amplitudes  in  aa  to  be 
larger  in  even  SC  DN  phases  compared  with  AX  or  odd  SC 
DN  phases.  We  found  Vsw  13.5-day  amplitudes  were 
higher  in  DN20  compared  with  DN21,  but  the  other  solar 
rotation  amplitudes  for  Vsw  and  Pe  were  generally  large  in 
all  DN  phases.  In  summary,  the  magnitudes  and  periodi¬ 
cities  of  Vsw  3nd  Pe  were  generally  largest  in  all  D  phases 
where  HSS  dominate.  The  expected  variations  in  the  Hale 
cycle  are  minimal,  in  part  because  the  22-23  Hale  cycle 
was  unusual. 

The  very  different  features  of  Vsw  and  Pe  in  every  solar 
cycle  are  thus  explained.  The  HSS  solar  velocity  sources 
tend  to  drive  the  observed  average  variations  of  both  the 
total  Vsw  and  the  total  Pe  and  their  solar  periodicities, 
although  they  only  account  for  approximately  half  of  the 
total  frequency  and  fraction  of  solar  structure  occurrences. 
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